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Abstract 



We consider the resummation of the logarithmic contributions to the region 
of small transverse momenta in the distributions of high-mass systems (lepton 
pairs, vector bosons, Higgs particles, . . . ) produced in hadron collisions. We 
point out that the resummation formulae that are usually used to compute 
the distributions in perturbative QCD involve process-dependent form factors 
and coefficient functions. We present a new universal form of the resummed 
distribution, in which the dependence on the process is embodied in a 
single perturbative factor. The new form simplifies the calculation of non- 
leading logarithms at higher perturbative orders. It can also be useful 
to systematically implement process-independent non-perturbative effects in 
transverse-momentum distributions. We also comment on the dependence of 
these distributions on the factorization and renormalization scales. 
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The properties of the transverse-momentum distributions of systems of high mass Q 
produced at high-energy hadron coUiders are important for QCD studies and for physics 
studies beyond the Standard Model (see, e.g., Refs. On the theoretical side the 

computation of these distributions is complicated by the presence of large logarithmic 
contributions of the type \n.Q'^ /q]_, which spoil the convergence of fixed-order perturbative 
calculations in the region of small transverse momenta q±. The logarithmically-enhanced 
terms have to be evaluated at higher perturbative orders, and possibly resummed to all 
orders in the QCD coupling ag. The all-order resummation formalism was developed in 
the eighties [P-||12[. The structure of the resummed distribution is given in terms of a 
transverse-momentum form factor and of process-dependent contributions. 

The transverse-momentum form factor is sometimes supposed to be universal, that is, 
independent of the process. Performing a process-independent calculation at relative order 



Og, two of us have recently shown in an explicit way that the form factor depends on the 
process. In this note, we give a general interpretation of the results in Ref. We point 



out that this process-dependent features persist to higher orders, at least in the final form 
in which the resummed transverse-momentum distribution is usually organized. We also 
present a new form of the resummed distribution, where all the logarithmic contributions 
are embodied in universal (process-independent) factors. 

We consider the inclusive hard-scattering process 

hi{p,) + h2{p2)-F{Q\ql)+X , (1) 

where the triggered final-state system F is produced by the collision of the two hadrons hi 
and /i2 with momenta pi and p2i respectively. We denote by -y/i the centre-of-mass energy 
of the colliding hadrons (s = {pi + P2Y — '^PiP2)- The final state F is a generic system 
of non-strongly interacting particles [|, such as one or more vector bosons (7*, W, Z, . . .), 
Higgs particles (H) and so forth. 

To simplify the discussion we limit ourselves to the case in which only the total invariant 
mass Q and transverse momentum (with respect to the direction of the colliding hadrons) 
of the system F are measured. The extension to more general kinematic configurations, 
in which, for instance, the rapidities are measured, is straightforward (see the discussion 
below Eq. (pS])). We also assume that at the parton level the system F is produced with 
vanishing q± (i.e. with no accompanying final-state radiation) in the leading-order (LO) 
approximation, so that the corresponding cross section is da'^^'^^ / dq\ oc 5(g^). Since F is 
colourless, the LO partonic subprocess is either qq annihilation, as in the case of 7*, and 
Z production, or gg fusion, as in the case of the production of a Higgs boson H. 

The transverse-momentum cross section for the process in Eq. (|I|) can be written as 



da I 



dQ^ dq] 



da I 



dQ^ dqj 



+ 



dai 



dQ^ dq] 



(2) 

fin. 



Both terms on the right-hand side are obtained as convolutions of partonic cross sections 
and the parton distributions fa/hix^Q"^) (a = qf,qf,g is the parton label) of the colliding 

^We do not consider the production of strongly interacting particles (hadrons, jets, heavy quarks, ...), 
because in this case the resummation formalism of smal\-q± logarithms has not yet been fully developed. 
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hadrons]^ The distinction between the two terms is purely theoreticaL The partonic cross 
section that enters in the resummed part (the first term on the right-hand side) contains 
all the logarithmically-enhanced contributions ag In*" Q^/g^.- Thus, this part has to be 
evaluated by resumming the logarithmic terms to all orders in perturbation theory. On 
the contrary, the partonic cross section in the second term on the right-hand side is finite 
order-by-order in perturbation theory when g_L — > 0. It can thus be computed by truncating 
the perturbative expansion at a given fixed order in as- 

The finite component of the transverse-momentum cross section is obviously process- 
dependent, and we have nothing to add on it in this paper. In the following we discuss the 
structure of the resummed part. 

The resummed component is0 

1 rl 



dap 



dQ^ dq] 



/ dxi dx2 db - Joibqi_) fa/hi{xi,bl/b'^) fb/h2ix2,bl/b'^) 
, Jo Jo Jo ^ 



■ W^,{xiX2s;Q,b) . (3) 

The Bessel function J^ibq^) and the coefficient b^ = 2e~'^^ {'Je = 0.5772 ... is the Euler 
number) have a kinematical origin. To correctly take into account the kinematics constraint 
of transverse-momentum conservation, the resummation procedure has to be carried out in 
the impact-parameter space. The transverse-momentum cross section @ is then obtained 
by performing the inverse Fourier (Bessel) transformation with respect to the impact pa- 
rameter b. The factor W^j^ is the perturbative and process-dependent partonic cross section 
that embodies the all-order resummation of the large logarithms In Q'^b^ (the limit q± Q 
corresponds to Qb ^ 1, because b is the variable conjugate to qj_). 

The resummed partonic cross section is usually (see, e.g., the list of references in Sec- 
tions 5.1 and 5.3 of Ref. |T|) written in the following form: 

W[,{s; Q,b) = f dz, f dz2 C^Mbl/b'), z,) CgiaM/b'), Z2) 6{Q' - z,Z2s) 
^ Jo Jo 

■ 4'''^''ms!{Q,b) . (4) 

Here, dcg*^'*'^ is the cross section (integrated over q±) for the LO partonic subprocess c+c 
F, where c,c = q,q (the quark and the antiquark qf can possibly have different fiavours 
/, /') or c, c = g, g. The expression (y\^^^^ can include an overall factor ag(Q^), as in the 
case of g + g ^ H through a triangular quark loop where p = 2. The term {Q, b) is the 
quark (c = q) or gluon (c = g) Sudakov form factor. The resummation of the logarithmic 
contributions is achieved by exponentiation that is by showing P, |^ that the form 

factor can be expressed as 



SciQ, 6) = exp J - [ 

( Jbl 



df_ 



A,(as(g^)) ln^ + i?,(as(g2)) 



(5) 



■''Throughout the paper we always use parton densities as defined in the MS factorization scheme and 
as (5^) is the QCD running couphng in the MS renormalization scheme. 

^As discussed at the end of the paper, this expression can be generahzed to include the dependence on 
the renormalization and factorization scales /i/j and [ip. 
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with c = q OT g. The functions Ac{as), Bc{as), as well as the coefficient functions Cabic^s, z) 
in Eq. contain no InQ^^^ terms and are perturbatively computable according to the 
power expansions^] 

oo 

n=l 

oo 

S.(as) = E(v)"^'"' • 

n=l 

oo 

Cab{as,z) = 5.b5{l-z) + Y,[^Yc^:^{z) . (8) 



TT 

n=l 



The knowledge of the coefficients A^^^ leads to the resummation of the leading logarithmic 
(LL) contributions. Analogously, the coefficients {y4^^\ C^^^} give the next-to-leading 
logarithmic (NLL) terms, the coefficients {A'^'^\ B'^'^\C^'^^} give the next-to- next-to-leading 
logarithmic (NNLL) terms, and so forth. The coefficients A^^\ A^'^\ B^^^ are known both 
for the quark |^ and for the gluon form factors 



= Cf , A« = Ca , 
Af = \cpK , Af = \caK , (9) 



where 



V18 6 y 9 ^ ^ ' 

The best studied example among the processes in Eq. (|I]) is lepton-pair Drell-Yan (DY) 
production through the LO partonic subprocess g + g ^ 7*(or W, Z) — * Z + V . In this case 
the NNLL coefficient E^^) was computed by Davies and Stirling [O]: 



^ 16 7 V36 48 2^7 ^ 18; 



where Cn is the Riemann (^-function (Q^ = 1.202 . . . ). The coefficient i?^^-* for Higgs boson 
production has recently been computed by two of us [|13]. The corresponding LO partonic 
subprocess is gluon fusion, g+g — > H, through a massive-quark loop. In the limit of infinite 
quark mass the value of -B^^-* is |1^ 



j^(2)H ^cl(- + -n'^- -Cs] + -CpNf -CANf(- + -]~ -CfCa . (12) 
^ ^ V24 18 2^7 2 ^ ■^V12 9y 8 ^ ' 

The main issue that we want to discuss in this paper regards the process dependence 
of the various factors in the resummation formula (^. As denoted by the superscripts F 

^Note that in Refs. ||l2|, ^ the perturbative coefficients are normahzed to powers of as /in rather than 
as/n. 
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in the various terms on the right-hand side, the coefficient functions depend on the 
process. This is conffimed by the calculations of the coefficients C^^h^ ^ performed in the 



literature for several processes [] [|T^] The form factor SdQ, b) that enters Eq, 



is (often) supposed to be universal (this is the reason why it is named quark or gluon form 
factor rather than DY, 77, WZ, H, . . . form factor). However, this is not the case: the 
form factor S^{Q,b) in Eq. (^ is process- dependent. In the following, we ffist present a 
universal (process-independent) version of the resummation formula (^ and we sketch its 
physical origin. We then clarify the relation between Eq. (P and our process-independent 
version. 

The process-independent resummation formula is 

W^,{s-Q,b) = 5^ / rf^i / dz2C,a{as{bl/b^),zi)C,b{as{bl/b^),Z2)6{Q'-ziZ2s) 
^ Jo Jo 

■ a^,{Q\as{Q')) S,{Q,b) . (13) 



It formally differs from Eq. (|) by the replacement (t[^^'' ^ {Q"^) — » (t^((5^, as(<5^))- While 
c"c5*^''^(Q^) is the cross section for the LO partonic subprocess, (7^(Q^, as(Q^)) includes 
higher-order QCD corrections to it, according to 

af,(Q^as(Q')) = ^''(g^) H!{as{Q')) , (14) 

where the function H^{as) has a perturbative expansion similar to Eqs. (§1)-(|): 

00 

J?f(as) = l + E(v)"*'""' • (15) 



TT 

n=l 



Note that the function (as) depends on the process. Nonetheless, its introduction is 
sufficient to transform the process-dependent form factor and coefficient functions C^^ 
of Eq. (|) into the process-independent form factor Sc and coefficient functions Cca of 
Eq. (PD. 



The resummation formula in Eq. (0), which can be derived by the customary resum- 
mation methods [|ll[], has a simple physical origin. When the final-state system F is 



kinematically constrained to have a small transverse momentum, the emission of accompa- 
nying radiation is strongly inhibited, so that only soft and collinear partons (i.e. partons 
with low transverse momenta qt) can be radiated in the final state (Fig. |I]). The process- 
dependent factor (as(Q^)) embodies hard contributions produced by virtual corrections 
at transverse-momentum scales qt ~ Q. The form factor Sc{Q, b) contains real and virtual 
contributions due to soft (the function ^^(as) in Eq. (^) and flavour-conserving collinear 
(the function Bc{as) in Eq. (§)) radiation at scales Q^qt^l/b. At very low momentum 
scales, qt^l/b, real and virtual soft-gluon corrections cancel because the cross section is 
infrared safe, and only real and virtual contributions due to collinear radiation remain 
(the coefficient functions Cab{o:s{bl/b'^), z)). Note that Sc{Q,b) and Cabiasibl/b"^), z) are 
process-independent and only depend on the flavour and colour charges of the QCD partons. 

"a general expression for the coefRcients C^V ^ in terms of the one- loop matrix element of the corre- 
sponding process is given in Eq. (17) of Ref. |13|. 
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Figure 1: Diagrammatic structure of the various factors that enter the process-independent 
resummation formula j^TBj). 



The two versions (^) and (0) of the resummation formula can formally be related as 
follows. We use the renormalization-group identity 



g,{as{Q')) = exp | T ^ ^72(«s(g')) 1 giiMbl/b')) 



(16) 



which is valid when 



92{as) = P{as)—-. , 17 

a In as 



where P{as) is the QCD /^-function 

din as{q'^ 



2 -f3{as{q')) , 



/3(«s) = -/3o— ) +... , 12Po = llCA-2Nf . (19) 

Then, setting gi^asiQ"^)) = H^^asiQ"^)) inserting the right-hand side of Eq. (|16D in 
Eq. (|l^), we immediately obtain Eq. (^. More precisely, the process-independent resumma- 
tion formula in Eq. (|T3|) implies the customary version in Eq. (^, provided the perturbative 
function Bc{as) in the form factor (see Eq. (|)) and the coefficient functions are 
related to their process-independent analogues by the following all-order relations 

CSXas.z) = [fff(as)]"' C„j(ag,j) , (20) 

BfM = B.K)-/3M '""f"'°^' . (21) 

dm as 

While the perturbative function Ac{as) and the first-order coefficient Bc^^ of the function 
-Bc(tts) are process-independent, the result in Eqs. ( pT]) and (|20|) shows that the coefficients 
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dependence of the first few coefficients is explicitly given by 



Be''' , Bc^' , . . . and C^^ (z), (z), . . . in Eq. (|) do depend on the process. The process 



c. 



(22) 

(23) 
(24) 
(25) 



Note that the process dependence of the resummation formula (^ is not simply embod- 
ied in pure numerical coefficients. For example, considering the Higgs production mecha- 
nism g + g ^ H through a massive-quark loop, the corresponding coefficient Hg and, 
hence, the form-factor coefficient Bg depend on the mass of the quark in the triangular 
loop. Moreover, the results discussed so far can straightforwardly be extended to more 
general configurations, in which several kinematical variables (and not only the invariant 
mass Q) of the final-state system F are measured. For example, when F is a pair of vector 
bosons, we can consider the corresponding transverse-momentum cross section at fixed ra- 
pidities of the vector bosons. The extension simply amounts to including the dependence 
on the kinematics of the system F in the cross section factor (and, hence, in cr^^^ and 
iff^) of the pro cess- independent resummation formula (|13|). In this general case, Eqs. (pl|) 
and (pOD imply that the form factor and the coefficient functions in Eq. (^) depend 
on the kinematics of the system F in a non-trivial manner. 

The process-independent resummation formula ([13|) can alternatively be used to relate 
the factors on the right-hand side of Eq. (^) for different processes. Considering Eqs. (pl|) 
and (|20|) for two different processes F and F, we obtain 



C|!(as,2;) 
5f(«s) 



^r(«s) 



1/2 



5,^(as)-/5(«s) 



Cc^("s,2;) , 
rflni7f^(as) 



din as 



(26) 
(27) 



where H^^ = /H^ . The perturbative expansion of these equations obviously leads to 
relations that are similar to those in Eqs. 



a 



cb \ 



Cir(.) + 5.,5(l-.)^i7«^^ 



^(2)F^^(2)F^^^^(1)FF ^ 

5(3)^ = i?(3)^ + /3ii/«^^ + 2/3o 



jji2) FF f ^(1) FF^ ^ 



(28) 

,(29) 

(30) 
(31) 



Since the function H^^ does not appear in Eq. (HI), within the framework of the process- 



dependent resummation formulae, Eqs. (^6D -(pl|) have to be regarded as non-linear relations 
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between the form factors and coefficient functions of different processes. For instance, 
having computed the first-order coefficient C^^'* for two processes, one can check Eq. 

and extract h'^^^^ . Then, the second-order coefficients B^^ of these processes must be 

related by Eq. (|30D . This constraint can either be checked or used to compute -BF'*^ from 



The relation in Eq. ( pS]) (or Eq. (0)) can straightforwardly be checked by comparing the 
known first-order coefficients C'^^b^ ^'^) production of DY lepton pairs [|1^, p!7[] , 

diphotons and ZZ pairs [^. In particular, we have 



^WFDY ^ jjiDF _ ^ii)DY = ]^Cf[V^ -7i^ + ^) , F = 77, ZZ , (32) 



where V"''^ and V^^ are given in Eq. (11) of Ref. []T8[ and in Eq. (8)|^ of Ref. [0, respectively. 



The relation (or (|D) for the NNLL coefficient B^^^^ has explicitly been derived by 
the process-independent calculation at Oi^a^) performed in Ref. [0. 

The reasoning used to obtain Eq. (H) from Eq. (|13D can also be used to show that 
Eq. (|r^ is invariant under the transformation 

H!(as(Q')) - fff(Qs«')) [j(as(Q'))]"' , 

SAQ,b) SJ,Q,b) expj r %B(a,(e)) "'""'"ff" ) . (33) 

\jbl/b^ r dlnas(r) J 

Ca,{as{hllh'),z) ^ a.(«s(&^/&'),^) b(«s(&^/&'))]'^' , 



where g{as) = l + 0{as) is an arbitrary perturbative function. This renormalization-group 
symmetry of the resummation formula ([13|) implies that its factors, H^, Sc (more precisely, 
the function Be) and Cab, although they are process-independent, are not unambiguously 
computable (defined) order by order in perturbation theory. The same conclusion can be 
reached by expanding Eq. (|13]) in powers of ag. At any given order in as, the coefficients 



of the various powers 

of lnQ252 depend on the unknowns i/i")"", Bi"\ but the number 
of such coefficients is less that the number of unknowns. For instance, at the first relative 
order in as, the coefficient of In^ Q^b'^ determines A^^\ the coefficient of InQ^b^ determines 
B^^\ while the coefficient of the constant term is a linear combination of H^^^^ and C^f^ . 

The perturbative ambiguity of the decomposition of the right-hand side of Eq. (0) in 
the factors H^, Sc and Cab is a consequence of the fact that the transverse-momentum cross 
section is not a collinear-safe quantity. The effect of coUinear radiation at low transverse- 
momentum scales, qt^l/b, is divergent in perturbation theory. The (arbitrary) regular- 
ization procedure of these divergences introduces some ambiguity in the definition of the 
coefficient functions Cab- Then, the ambiguity propagates to Sc and through collinear 
evolution (see Eq. (0)). In this way of thinking, the ambiguity is similar to that encoun- 
tered in the definition of the parton densities. As the parton densities have to be defined 
by fixing a factorization scheme (e.g. the MS scheme or the DIS scheme), the factors on the 
right-hand side of Eq. (|1^) has to be defined by choosing a 'resummation scheme'. Note 
that the choice of a 'resummation scheme' amounts to defining (or Cab) for a single 



*As pointed out in Ref. [13|, the actual expression in Eq. (8) of Ref. ||l£[| is not correct. 
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processQ. Having done that, the pro cess- dependent factor and the universal factors Sc 
and Cab in Eq- ( P^ are unambiguously determined for any other process in Eq. (0). 

We can suggest three examples of 'resummation schemes', that is three prescriptions 
to use the process-independent factorization formula ([T3| ) in practice. A first 'short-cut' 
(because the DY transverse-momentum distribution is best studied) prescription, which we 
name 'DY resummation scheme', is obtained by setting if^^(as) = 1 in the case of the DY 
cross section (integrated over the rapidity of the vector boson). This unambiguously fixes 
the process-independent form factor Sc and coefficients functions Cab as those determined 

(1') 

from the DY process. In particular, the coefficient Bq of the quark form factor is that 
in Eq. (|ll|), and the quark coefficient functions C^^^^[z) {b = q,q,g) are those computed 
in Refs. 14, |l5|, for the DY process. Any other gg-initiated cross section is then 



obtained by using Eq. (|1^) with a perturbatively computable hard function H^{as)- This 
prescription can be extended to ^f^f-initiated cross sections by setting H^{as) = 1 in the 
case of Higgs production via gg fusion. Thus, in the case of an infinitely-massive quark in 
the loop (the definition of the scheme will be different by keeping the mass of the quark 
finite), the coefficient Bg^^ of the gluon form factor is that given in Eq. (|12D and the gluon 
coefficient functions C'^]^{z) {b = q,q,g) are those computed in Ref. [|l^ . 

A second prescription, which is in some sense more physical, exploits the fact that the 
first moment (with respect to Q'^/s) of the flavour non-singlet (NS) hadronic (and partonic) 
cross sections deflnes collinear-safe quantities. This property follows from fermion-number 
conservation]]]. Thus, we can deflne a 'NS resummation scheme' by flxing the overall nor- 
malization of the NS coefficient function Cab{(Xs,z) in Eq. (|1^) in such a way that its first 
moment (with respect to z) vanishes. In this scheme (the subscript NS labels the scheme 
choice) 

1 

dz [ Cgq^Ns{as, z) - Cgq^Ns{as, z) ] = 1 , (34) 



and thus, using the DY coefficients functions from Refs. |T^, |T^, |1^, and Eq. (0), we 
obtain 

rii)DY ^ /ti"^ 7 



Kns = ^ Y - 2 J • ^35) 

Inserting Eq. (|35|) in Eq. (^4]), and using the DY coefficient in Eq. ([TT|), we then obtain the 
corresponding coefficient of the quark form factor: 

r(2) _ r>(2) DY n tt{1)DY 

The reason why this scheme can be considered 'more physical' is that the first moments 
of the NS cross sections are collinear-safe and, hence, free from ambiguities related to the 



*More precisely, H^^ has to be defined for two processes: one process tliat is controlled, at LO, by qq 
annihilation and another process that is controlled, at LO, by gg fusion. 

Tin the case of the Altarelli-Parisi evolution of the parton densities, fermion-number conservation implies 
the vanishing of the first moment of the Altarelli-Parisi probabilities. 
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regularization procedure of coUinear singularities. As is known 11, 20 1, the process- 
independent form factor coefficients in Eq. @ measure the physical intensity of soft {A^c^ 
and A^c^) and collinear (-Bc^^) radiation. Since the coefficient B^^j^g in Eq. (|3^) appears in 
the form factor of an infrared and collinear safe quantity, it should measure the intensity 
of the collinear radiation from quarks at the second order in as and it should enter in 
the form factor of other infrared and collinear safe quantities. For example, an analogous 
transverse-momentum form factor controls the back-to-back region of the energy-energy 
correlation (EEC) in e~^e~ annihilation ^T|. The coefficient in Eq. ( p6D should directly be 
related to the corresponding coefficient B^'^^ of the EEC, modulo (possible) effects due to 
the difference between the space-like kinematics of the transverse-momentum cross sections 
and the time-like kinematics of e^e~ annihilation. 

The generalization of the 'NS resummation scheme' to ^f^f-fusion processes can be ob- 
tained, for instance, by fixing the overall normalization of the gluon coefficient function 
Cgg{as, z) in Eq. (pTSD in such a way that its first moment vanishes. Note, however, that in 
the gluon (or, more generally, flavour- singlet) channel there is no anologue of the fermion- 
number conservation rule. Flavour-singlet transverse-momentum cross sections are collinear 
unsafe quantities. Thus, a 'physical' (collinear-safe) interpretation of the higher-order co- 
efficients of the function Bg{as) of the gluon form factor is less straightforward than in the 
quark channel. 

The coefficient B^'^ in Eq. (^) coincides with the coefficient of the end-point contribution 
(the term proportional to 5(1 — 2;)) to the LO Altarelli-Parisi probability Pcc\z). As can 
be noticed from Eqs. (22) and (23) in Ref. |13|, the coefficient B^J^g in Eq. (^) (and the 
coefficients Bjj^^^^ and Bg^^^ in Eqs. (|TT]) and (|T2|)) does not coincide with that of the end- 

point contribution to the NLO Altarelli-Parisi probability Pec (z) in the MS factorization 
scheme. This should not be surprising, because beyond LO the Altarelli-Parisi probabil- 
ities are not collinear-safe quantities: they depend on the regularization procedure of the 
collinear singularities and on the factorization scheme. This interpretation is confirmed by 
the fact that the various coefficients Be differ from one another by terms proportional to 
the first coefficient /Sq of the /3-function, as should be expected from factorization-scheme 
dependence. Using Eqs. (|2D|)-(|2S|), it is obviously possible to introduce an 'MS resum- 
mation scheme' such that the function Bc{as) of the process-independent form factor in 
Eqs. (^ and (|T3|) coincides by definition with the perturbative function B^jj^{as) that 



controls the end-point contribution to the all-order Altarelli-Parisi probabilities in the MS 
factorization scheme. 



Note that our discussion on the relation between Eqs. (^) and (|T3D and on the resummation- 
scheme ambiguity does not imply that the two equations are equivalent. In fact, the 
process-independent resummation framework contains more information than its custom- 
ary process-dependent version. To evaluate resummed cross section within the framework 
of Eq. (^), we should compute two functions (S'f and C^) for each process F. Once a 
'resummation scheme' has been defined, a similar evaluation by using Eq. (|I3|) requires the 
computation of two universal functions {Sc and Cab) and of a single additional function 
(H^) for each process. 

Of course, we can still continue to use Eq. M). In this case, we can exploit the addi- 
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tional information contained in Eq. (13) by using Eqs. (26) and (27) (or their perturbative 
expansions in Eqs. (p^)-(|3T|)) to relate form factors and coefficient functions for different 
processes. 

Independently of the resummation formula that is actually used, once the form factor 

(2) 

coefficient Be has been computed for a single process by performing a calculation at 
relative order ag, no further (9(ag)-calculation is necessary for its universal implementation 
in all processes. The implementation can simply be performed by computing the coefficient 
functions at 0{as) and then using Eqs. (^) and ( plf ) (or Eqs. (^) and (|30D). Analogously, 
using Eqs. (^) and ( ^5]) (or Eqs. (^) and (0)), the universal implementation of the form 
factor coefficient si^^ requires its 0{al) computation in a single process and the evaluation 
of the coefficient functions at 0{a^) for the various processes. This procedure extends to 
any higher orders. 

At present, owing to the results in Eqs. (|Il|) and (|12|), the NNLL coefficient Be can 
be included in all calculations of transverse-momentum cross sections in hadron collisions. 
Its inclusion cannot be regarded as the extension of these resummed calculations to full 

(3) 

NNLL accuracy, because the corresponding coefficient Ac is not yet known. Nonetheless, 

(2) 

the knowledge of Be can certainly be used to improve the matching with fixed-order 
calculations at high q±. The matching can be performed by supplementing the resummed 
cross section with the finite contribution on the right-hand side of Eq. (^. The finite 
contribution is obtained, from the complete calculation at a fixed order in ag, by subtraction 
of the corresponding perturbative terms already included in the resummed component. If 
B^^^ is included in this component, the remaining finite component at relative order ag 
contains no InQ^/g^. terms. Thus, the latter can be inserted in Eq. (|^) uniformly with 
respect to q±, that is without including logarithmic terms that spoil the convergence of the 
matched result in the small-g^ region. 



In this paper we have not discussed how non-perturbative effects (see Refs. [^-[^ and 
references therein) affect the transverse-momentum cross sections. Although the process- 
independent resummation formula (p!3D can be recast in the form of Eq. (^ at the pertur- 
bative level, it is not evident whether the two formulations are equally suitable to deal with 
non-perturbative contributions. In particular, since the form factor and the coefficient 
functions in Eq. (H) are process-dependentj^, a simple and uniform implementation of 
process-independent non-perturbative effects (such as those due to the initial intrinsic kt 
of the partons in the colliding hadrons) in Eq. (^) can be very involved. The process- 
independent version (^) of the resummed cross section can help to consistently introduce 
and implement non-perturbative effects in transverse-momentum distributions. 

We conclude by discussing the dependence of the resummation formulae on the renor- 
malization and factorization scales and fip. This dependence is often parametrized by 
introducing some arbitrary coefficients, Ci,C2,C3, as suggested in Refs. [T^. This is 
a perfectly sensible and reasonable procedure to try to estimate the effect of higher-order 
corrections, but it does not exactly correspond to the procedure that is usually followed in 
fixed-order perturbative calculations. In the case of soft-gluon resummed calculations for 

^We remind the reader that the process dependence can kinematically be quite non-trivial when the 
transverse-momentum cross section in not fuUy integrated over the kinematical variables of the final-state 
system F. 
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event shapes and jet rates in e+e" annihilation and for threshold contributions in hadronic 
cross sections, a different procedure was introduced in Refs. [^] and . The latter proce- 
dure closely matches the use of and fip fixed-order calculations, and it can directly be 
applied (see, e.g., Refs. ^ and several other references quoted in Section 5 of Ref. 
to compare the scale dependence of resummed and fixed-order predictions. Here, we would 



like to note that the procedure of Refs. [l28| , |29[| can be introduced also in the case of the 
transverse- momentum distribution^ 

To this purpose, we first perform the replacement fa/h{x, h'^/lP') fa/h{x, in Eq. (|^) 
and we rewrite the resummed component of the transverse-momentum cross section as 



dap 



oo ^ 



To obtain Eq. 
given by 



^ / dXi I dX2 fa/hA^l, fb/h2{x2, 

W[i,{xiX2s;Q,b,ij,F) ■ (37) 
we have simply used the scale dependence of the parton densities as 

^ dz 



db - Jo( 



1 dz 

fa/h{x,bl/b'^) = ^ — Uabiz;bl/b^, fXp) fb/hix/z, f/p) , 



(St 



where Uab{z; b^/b"^, fij^) is the customary evolution operator matrix obtained by solving the 
Altarelli-Parisi evolution equations to the required logarithmic accuracy. The resummed 
partonic cross section in Eq. (|37| ) is thus obtained from that in Eq. (|l^) (or (^)) by re- 
placing the coefficient functions Cca{ois{bl/b'^),x) (or C^^{as{bl/b'^),x)) with the following 
convolution 

/•I ^7 

CUcysibl/b'),x;fil) = J2 -C,b{cys{bl/b'),z)UUz;bl/b^f^l) . (39) 

, Jx Z 



Note that, unhke the parton densities in Eq. (^, those in Eq. ( pT] ) do not depend on the 
impact parameter b, but only on the factorization scale fip. Expression (^) can thus be use- 
ful in practical calculations, because it avoids the numerical integration of the parton den- 
sities with respect to the impact parameter. The only factor of the resummed cross section 
that has to be integrated over b is now the perturbative component W^{xiX2s; Q, b, fip)- 



+ ... 



Then, to proceed as in Refs. |28, 2^, we simply observe that we can choose fip and 
fiR to be of the order of the hard scale Q. Thus, W^f^{xiX2s; Q,b, ftp) contains only an 
additional large parameter, L = InQ^fe^, and it can be written as 

W^{Q,b,HF) ~ exp (l /i fas^) + /2 ( "s^; — , — ) + «s/3 ( "s^; — , — 

■ a^(Q^as;QV/i?^) , (40) 

where as = as(yu|.). The notation on the right-hand side is symbolical, because we have 
understood the dependence on the momentum fractions [] 2:1, and on the parton indices 

§A similar observation can be found in Ref. [2^ . 

^The convolutions with respect to the momentum fractions can as usual be diagonalized by taking 
A^-moments and working in Mellin space. 
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a,b = q, q, g. In particular, the exponential should be understood as an exponential matrix 
that depends on the parton indices. 

The form in Eq. (|40|) is straightforwardly obtained by using the exponentiated result in 
Eq. for the form factor and the customary exponential form for the evolution operator 
Uab{z; bl/b"^, fj,p). In the exponent, the function Lfi resums all the LL contributions a^L^~^^, 
the function /2 resums the NLL terms agL", while as/s contains the NNLL terms a^L^~^, 
and so forth. The functions are normalized as fi{asL = 0) = and can easily be obtained 
in analytic form. 

Although it is symbolical, the form of Eq. (^0|) is sufficient to make our point. Indeed, 
it is exactly in the same form as the expressions that are obtained by performing soft-gluon 
resummation for e~^e~ event shapes and hadronic cross sections near threshold. Recasting 
the resummed transverse- momentum distributions in the form of Eq. (^OD, the theoretical 
accuracy of the resummed calculation can be investigated as in fixed-order calculations, 
by varying and fiji around the value Q of the typical hard scale. Future studies along 
these lines (for instance, comparisons of the scale dependence of resummed and fixed-order 
calculations) could be useful to increase our confidence in the quantitative rehabihty of 
resummed predictions for transverse-momentum cross sections. 
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